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ABSTRACT 


In  the  past  year  I  have  been  processing  a  large  amount  of  regional/teleseismic  data  from  various  broad¬ 
band  seismic  stations  in  eastern  Eurasia.  Eourier  spectra  of  Pn,  Eg  and  Pg  waves  were  computed  for 
many  events  and  paths  to  study  path  attenuations.  Among  the  data  collected  and  processed  are  Eg 
spectra  collected  from  the  1991-1992  Tibetan  Plateau  Passive  Experiment.  Using  these  spectra  and  a 
standard  two-station  method  that  virtually  eliminates  source  and  site  effects,  I  obtain  a  model  of 
Qo  =  (126  ±  9)  and  j]  =  (0.  37  -I-  -0. 02)  in  a  frequency  range  between  0.2  and  3.6  Hz,  where  Qq  and  rj  are 
Eg  Q  at  1  Hz  and  its  power-law  frequency  dependence,  respectively.  The  estimated  Qq  value  is  among  the 
lowest  ever  reported  for  continental  areas;  it  qualitatively  supports  the  observation  by  McNamara  et  al. 
(1996)  that  Eg  can  not  be  observed  inside  the  plateau  beyond  about  700  km,  a  limiting  distance  that  is  much 
shorter  than  those  in  the  other  low  Qq  (~  200)  regions,  such  as  Iran  and  the  western  U.S. 

Quantitatively,  the  estimated  Qq  value  is  lower  by  a  factor  of  3  than  the  values  of  366  estimated  by  McNa¬ 
mara  et  al.  (1996),  who  used  data  from  the  same  experiment.  Since  there  are  several  differences  in  the  data 
processing  and  inversion  procedures  used  in  this  and  the  previous  studies,  I  investigated  the  effects  of  these 
differences  on  the  Q  estimates.  I  conclude  that  the  most  probable  cause  of  the  discrepancy  is  in  the  differ¬ 
ent  inverse  methods  used.  This  is  so  because  the  previous  inversion  allowed  the  source  and  site  terms  to  be 
free  parameters  solved  for.  Since  20  events,  8  stations  and  5  frequency  bands  were  used,  the  unknown 
source  and  station  terms  should  be  more  than  100.  In  this  study  only  two  free  parameters  (Qq  and  rj)  are 
solved  for,  thus  avoiding  the  instability  caused  by  parameter  trade-offs. 

This  research  suggests  that  the  previously  observed  Eg  blockage  for  paths  crossing  the  northern  boundary 
of  the  plateau  may  be  partially  or  entirely  caused  by  the  abnormal  low  Eg  Q.  Further  research  on  lateral 
and  depth  variations  of  crustal  Q  in  and  around  the  Tibetan  Plateau  is  highly  recommended. 

KEY  WORDS: Eg  Q,  Tibetan  Plateau,  Eg  blockage.  Central  Asia. 


OBJECTIVE 


The  primary  objective  of  this  research  is  to  quantify  path  attenuation  of  regional  waves  in  continental  areas, 
such  as  Eurasia,  by  developing  digital  tomographic  Q  maps.  The  proposed  research  include  (a)  various 
improvements  of  the  methods  for  measuring  path-variable  Qq  and  ;/  (Q  at  1  Hz  and  its  power-law  frequency 
dependence,  respectively)  using  regional  wave  spectra,  (b)  applications  the  improved  methods  to  regional 
waves  paths  in  Eurasia  to  measure  Qq,  rj,  and  (c)  to  input  the  measured  Q  values  to  a  computerized  tomo¬ 
graphic  algorithm,  to  obtain  laterally  varying  Q  maps  for  Eg  and  other  regional  waves. 

This  research  provides  important  input  to  the  world-wide  monitoring  of  nuclear  explosions.  The  tomo¬ 
graphic  Q  maps  can  be  used  for  the  calculation  of  source  spectral  characteristics  of  any  future  seismic  event 
to  infer  the  nature  and  size  of  the  event.  These  Q  maps  can  also  be  used  for  estimating  the  detection 
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threshold  of  the  existing  seismic  stations  located  within  the  area  studied. 

RESEARCH  ACCOMPLISHED 

Over  the  past  I  have  continued  data  collection  and  inversion  for  Lg  Q  in  Eurasia.  In  this  paper  I  report  anal¬ 
ysis  of  a  data  set  that  consists  of  Lg  spectra  from  1 1  broadband  PASSCAL  seismic  stations  deployed  during 
the  1991-1992,  passive  Tibetan  Plateau  experiment  (McNamara  et  al.,  1996).  In  an  effort  to  invert  for 
source  spectra,  I  found  that  the  data  required  that  the  Lg  Q  in  eastern  Tibet  be  re-measured. 

A  robust,  two-station  method  is  used.  Ideally  the  method  uses  stacked  spectral  ratios  (SSR,  see  equations 
(10)-(17)  of  Xie  and  Mitchell,  1990)  that  are  calculated  using  two-station  pairs  that  are  aligned  to  the  same 
event-to-station  azimuths.  In  practice  the  SSRs  are  obtained  by  requiring  that,  the  differences  between  the 
azimuths  from  the  same  events  to  the  two  stations  are  smaller  than  a  preset  maximum  allowable  value, 
(‘5^)max-  The  choice  of  (<y0)niax  is  less  restrictive  for  the  Lg  than  for  many  other  phases,  since  the  Lg  con¬ 
tains  a  minimal  source  radiation  pattern  in  3D  structures.  In  this  study  two  different  values  of  of 

30°  and  12°  are  used.  These  values  result  in  37  and  22  two-station  pairs  from  the  data  set,  respectively. 
SSRs  are  formed  and  used  to  fit  Qigif)  models.  Ligure  1  shows  the  path  coverage  by  the  two  sets  of  SSRs, 
and  Ligure  2  shows  the  SSRs  and  the  fit  of  the  Qigif)  models.  SSRs  obtained  with  =  30  °  are  fit¬ 

ted  by 


Qigif)  =  (126  ±  9)/(0-37±o.o2)  q  2Hz<f<  3. 6Hz  . 

SSRs  obtained  with  a  of  12°  are  very  similar  to  those  with  a  (56)^^  of  30°  (Ligure  2),  and  are  fit¬ 

ted  by  a  Qigif)  with  Qq  of  (134  ±  10)  and  rj  of  (0.  32  ±0.02),  respectively.  These  values  are  virtually  the 
same  as  those  in  equation  (1),  but  are  theoretically  more  subject  to  random  errors  and  bias  owing  to  a 
smaller  dataset.  Qigif)  in  equation  (1)  is  therefore  my  preferred  model.  This  model  is  in  agreement  with 
the  short  limiting  distance  of  about  700  km  for  Lg  observation  in  the  study  area  reported  by  McNamara  et 
al  (1996),  and  the  recent  Qq  estimates  by  Lan  and  Lay  (2001). 

Xie  (2001)  discussed  in  detail  about  several  possible  reasons  why  the  Q  model  in  this  paper  differs  from  a 
previous  one,  and  concluded  that  the  most  likely  cause  is  a  parameter  trade-off  in  certain  inverse  schemes. 

CONCLUSIONS  AND  RECOMMENDATIONS 

The  Qq  value  of  about  126,  estimated  for  the  eastern  Tibetan  plateau  in  this  study,  is  among  the  lowest  ever 
documented  for  any  continental  areas.  An  implication  of  this  low  value  is  that  the  well-known  "blockage" 
of  Lg  for  paths  crossing  the  boundaries  of  the  plateau  (e.g.,  Ruzaikin  et  al.,  1977;  Ni  &  Barazangi,  1983) 
may  be  largely,  or  even  entirely,  attributed  to  the  low  Qq  values  in  the  plateau  (see  Ligure  3  for  a  detailed 
discussion).  Another  implication  of  the  low  Qq  is  that  the  crust  in  Tibet  may  be  characterized  by  higher- 
than-normal  temperature  and  fluid  content,  which  are  responsible  for  the  low  Qq  values.  These  are  in  line 
with  the  electral  conductivity  measurements,  and  suggests  that  Lg  Q  should  vary  laterally,  sometimes  dras¬ 
tically,  in  the  Tibetan  plateau. 

Luture  research  should  be  directed  to  analyzing  more  seismic  data  from  recent  seismic  experiments  in  the 
plateau,  to  resolve  details  of  the  lateral  variations  of  Qigif)  in  the  plateau. 
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Figure  1.  Locations  of  the  PASSCAL  stations  deployed  during  the  1991-1992,  passive  Tibetan  Plateau  experiment 
(solid  triangles),  earthquakes  (open  circles)  and  explosion  (star)  used  in  this  study.  Solid  paths  are  those 
satisfying  a  {56)^^^  (the  maximum  allowable  difference  between  the  event-to-station  azimuths  of  two  sta¬ 
tions;  see  text)  of  12°  when  two-station  pairs  are  selected.  Dashed  paths  are  those  satisfying  a  (<J6')max  of 
30°.  More  information  of  the  stations  and  events  can  be  found  in  McNamara  et  al.  (1996). 
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f  =  0.2-3.6hz,Qo  =  126  +  -9,ri=0.37  +  -0.02 


frequency(hz) 


Figure  2.  Stacked  spectral  ratios  (SSRs)  from  many  two-station  pairs  plotted  in  Figure  1,  and  the  fit  of  best  Q  mod¬ 
els  (straight  lines).  Black  and  gray  symbols  are  SSRs  obtained  using  a  of  30°  and  12°,  respec¬ 

tively.  The  Q  models  from  fitting  both  sets  of  SSRs  are  similar.  The  Q  model  written  on  the  top  of  the 
panel  is  from  fitting  the  black  symbols. 
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strong  scattering  at  the  topographic  boundary  is  not  required.  This  is  consistent  with  the  simulation  of,  e.g.,  Kennett  (1986),  who 
shows  that  the  scattering  is  unlikely  to  fully  account  for  the  blockage. 


